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1. Introduction 


While the distribution of neutral hydrogen mapped by SKA1-LOW provides an excellent and 
unique view of the reionization process over a large range of redshifts, detecting the sources re¬ 
sponsible for reionization directly sheds light on the crucial stage of galaxy formation and comple¬ 
ments our understanding of EoR. Extremely deep imaging with the Hubble Space Telescope (HST) 
has begun to probe the very bright end of the UV luminosity functions at z > 6 (Bouwens et al. 
2014; Robertson et al. 2013), with improvements expected with the James Webb Space Telescope 
(JWST). In the sub-nim, the Atacama Large Millimeter Array (ALMA) has detected individual 
high redshift, luminous objects known from existing surveys (e.g., Ouchi et al. (2013)). However, 
observations that are aimed at detecting individual galaxies at z > 6 arc difficult and time consum¬ 
ing, and neither of these space-borne facilities nor ALMA is expected to resolve the bulk of low 
luminosity sources responsible for reionization (Salvaterra et al. 2011). Approaches which can 
access the entire luminosity function of reionizing sources arc needed. 

Line Intensity Mapping has emerged as a promising technique that is sensitive to the integrated 
light produced by faint galaxies: instead of resolving individual sources, one measures on larger 
spatial scales the collective emission from an ensemble of sources, while retaining the spectral-thus 
redshift-information. This allows efficient redshift surveys that probe the integrated luminosity 
function of sources and provide three-dimensional information to study star formation activities 
during EoR. 

A few spectral lines arc currently being considered as promising tracers for high-redshift star 
formation activities in the intensity mapping regime. Among them, the most promising ones arc 
the rotational transitions from carbon monoxide (CO) (Righi et al. 2008; Visbal & Loeb 2010; 
Carilli 2011; Gong et al. 2011; Lidz et al. 2011; Pullen et al. 2013; Breysse et al. 2014), the 158jtm 
emission from singly ionized carbon ([CII]) (Gong et al. 2012; Uzgil et al. 2014; Silva et al. 2014), 
and the Lyman-a transition line from hydrogen (Silva et al. 2013; Pullen et al. 2014). Such large- 
scale surveys will not only reveal early star formation history but also measure the clustering of 
ionizing sources. These line intensity maps mark the three-dimensional distribution of ionized 
regions and probe different gas phases that complements the 21cm EoR surveys which trace neutral 
hydrogen. Together, they draw a complete view of the reionization process in the high-redshift 
Universe. In addition, on scales larger than the ionized regions, these line tracers anti-correlate 
with the 21cm emission. By measuring the shape of the cross-power spectrum of the two surveys, 
one can determine the characteristic scale of ionized regions, by marking the scale where the cross¬ 
correlation becomes negative, and be able to constrain statistically the characteristic size scale of 
ionized regions as a function of redshifts (Righi et al. 2008; Lidz et al. 2009). 

Below we discuss each of these tracers in detail, and present forecasts on the measurements 
of power spectra and cross-correlation signals with the SKA 1-LOW 21cm EoR survey. Due to 
the short emission wavelengths, it is not possible to observe [CII] and Lya with the SKA, thus 
we present the results assuming other future surveys. Lor CO, however, the proposed highest- 
frequency band of SKA 1-MID can potentially cover the redshifted CO(l-O) transition at z > 7.3. 
We discuss such possibilities in the next sections. 
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2. CO Intensity Mapping 

The CO(l-O) rotational line has a rest frequency of 115 GHz. The proposed highest fre¬ 
quency band of SKA1-MID, 4.6-13.8GHz, will have a chance of capturing the redshifted CO(l-O) 
at z > 7.3. For intensity mapping purposes, where we aim to measure the large-scale distribution 
of redshifted CO at low-angular resolution, only the inner few kilometer core of SKA1-MID will 
be relevant. The uncertainty in the theoretical modelling of CO brightness temperature at high 
redshifts is still large, but most predict the amplitude to be about 10 6 K or smaller on quasi-linear 
scales. Large collecting area and high sensitivity are required, and a densely packed antenna config¬ 
uration at the central core is desired, much like the requirement for SKA1-LOW for the 21cm EoR 
survey. Currently, however, the designed filling factor of SKA1-MID at the core is relatively small, 
only about 10 5 in the inner 1 km core in diameter, making the prospect of detecting redshifted CO 
very challenging. 

Here we instead use SKA 1-MID as a collection of single dishes, where each antenna records 
the total power while the cross products between antenna pairs are ignored. We note that perform¬ 
ing the CO survey in the single dish mode lifts the compact antenna configuration requirement 
and enables a survey in the intensity mapping regime, but imposes constraints on the stability of 
the auto-correlation spectrum which would not be present if it were done interferometrically. We 
consider a narrow redshift range at z = 8 ± 0.5, the lowest possible redshift allowed by SKA1-MID 
which should be more CO-rich. The specification is listed in Table 1, which also specifies the 
assumed SKA1-LOW 21cm survey parameters. Note that the assumed survey areas for CO auto¬ 
power spectrum and CO-HI cross-power spectrum calculations are chosen to be 0.1 and 10 deg 2 , 
respectively, in order to optimize the signal-to-noise ratio in each case. 

The amplitude of predicted CO brightness temperature fluctuations is rather uncertain, dif¬ 
fering by orders of magnitude at z = 8 from model to model. This is currently one of the most 
challenging aspects for CO intensity mapping work. Figure 1 shows the forecasted CO power spec¬ 
trum and the COx21cm cross-power spectrum measurements with the survey parameters listed in 
Table 1. Here we demonstrate the model uncertainties by plotting three of them: the middle solid 
curve with error bars is based on the models in Righi et al. (2008), while the upper dashed curve 
is based on Lidz et al. (2011) and lower dashed curve from Gong et al. (2011). The error bars 
include contributions from thermal noise and cosmic variance. The predicted total signal-to-noise 
ratio (SNR), summed over all accessible scales, of the upper, central and lower power spectra with 
SKA1-MID in single dish mode are 20.7, 4.0 and 0.1; for the CO and 21cm cross-power spectra 
with SKA1-MID and SKA1-LOW, the corresponding SNRs are 26.1, 6.0 and 0.7. The CO(l-O) 
signals and the cross power spectrum between 21cm and CO(l-O) emissions would be detectable 
at statistically significant level in more optimistic scenarios. The CO model uncertainties, how¬ 
ever, exceed the measurement errors, thus making it difficult to plan a survey. On the other hand, 
a detection of CO power spectrum can be very discriminating against models and guide our the¬ 
oretical understandings. Observational efforts for mapping CO at modest redshifts (z = 2 — 3) 
with existing instruments have only reached initial results, probing the high-end of CO luminosity 
functions (Decarli et al. 2014; Walter et al. 2014; Hodge et al. 2014). It is essential to improve 
our theoretical modelling and advance observational measurements of CO brightness temperature 
fluctuations across redshifts, which is currently an active area of research. In addition, one may 
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need to worry about foreground contaminations from synchrotron and free-free radiations coming 
from the Galaxy and other extragalactic sources. The severity depends on the observing frequen¬ 
cies (or the redshifted CO rotational lines of interest) and the CO signal strengths. However, these 
foregrounds are expected to be spectrally smooth and can be separated from the line emissions 
(Angelakis et al. 2009; Keating et al. 2014). At a few tens of gigahertz frequency range, spinning 
dust may also be a potential contamination, although its spectrum and strength are not well known 
(Ali-Halmoud 2013). It will likely require a dedicated experiment with both a larger field-of-view 
and greater surface brightness sensitivity to make a strong CO detection, or a significantly more 
compact SKA-MID core than currently envisaged. 



k (h/Mpc) k (h/Mpc) 


Figure 1: Left: SKA1-MID CO power spectrum at z = 8. Right: SKA1-LOW 21cm x SKA1-MID CO 
at z = 8. The upper, central, and lower curves indicate predicated CO signal strengths from three different 
models (see texts), whose uncertainties are larger than the predicted error bars. The survey areas for CO 
auto-power spectrum and CO-HI cross-power spectrum calculations are assumed to be 0.1 and 10 deg 2 , 
respectively, to optimize the signal-to-noise ratios. 


3. CII Intensity Mapping 

Carbon is one of the most abundant elements in the Universe and it becomes singly ionized 
[CII] with an ionization energy of 11.26 eV, less than that of hydrogen. With a splitting of the 
fine-structure level at 91 K, [CII] is easily excited resulting in a line emission at 157.7 fim through 
the ~Pi /2 —> 2 P \/2 transition. It is well established that the bulk of [CII] emission comes from 
photodissociation regions (PDRs), and provides a major cooling mechanism for the neutral inter¬ 
stellar medium (ISM). It is generally the brightest emission line in star-forming galaxy spectra and 
contributes to about 0.1 % to 1 % of the total far-infrared (FIR) luminosity in low redshift galaxies. 
Since carbon is naturally produced in stars, [CII] emission is expected to be a good tracer of the gas 
distribution in galaxies. ALMA high-resolution observations have revealed [CII] in high-redshift 
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Table 1: Parameters for SKA1-LOW and -MID at z=8±0.5. For SKA1-MID, single dish observation mode 
using 254 antennae is assumed. The survey areas for CO auto-power spectrum and CO-HI cross-power 
spectrum calculations are assumed to be 0.1 and 10 deg 2 , respectively, in order to optimize the signal-to- 
noise ratios. 



SKA1-FOW 

SKA1-MID 

auto/cross 

unit 

ant. diameter D ant 

35 

15 

m 

survey area A s 

13 

0.1/10 

deg 2 

FoV per ant. 

13 

0.01 

deg 2 

effective area per ant. A e 

925 

170 

m 2 

freq. resolution dv 

3.9 

9.7 

kHz 

bandwidth (z=8±0.5) BW 

18 

1427 

MHz 

tot. int. time t mt 

1000 

10,000 

hr 

min. baseline D m ; n 

30 

- 

m 

max. baseline D max 

1 

- 

km 

uv min 

16 

- 


uv max 

526 

- 


T 

A sys 

400 

25 

K 

effective num. ant. 

433 

254 


num. density of baselines «base 

0.8 

- 



galaxies, e.g, (Riechers et al. 2014; De Breuck et al. 2014), although no detections have been made 
for [CII] associated with galaxies at z >~ 7 (Gonzalez-Lopez et al. 2014; Ota et al. 2014) 

Even if the angular resolution to resolve the [CII] emission from individual galaxies at high 
redshift is not available, the brightness variations of the [CII] line intensity can be used to map the 
underlying distribution of galaxies and dark matter (Basu et al. 2004; Visbal & Loeb 2010; Gong 
et al. 2012). 

Here we follow Gong et al. (2012) to calculate the expected [CII] line fluctuations at 6 < z < 9, 
assuming that [CII] emission mainly originates in the hot gas in galaxies and that it is proportional 
to the gas metallicity, based on both analytical arguments and numerical models. Alternatively, 
[CII] emission can be estimated using observational relations, such as the empirical relation be¬ 
tween [CII] luminosity and star formation rate (SFR) in low redshift galaxies from Sargsyan et al. 
(2012). This has the advantage of including [CII] emission from several media, since [CII] is emit¬ 
ted not only from hot ionized gas but also from cold, mostly neutral gas in PDRs, and also in minor 
proportion from other regions. The connection between [CII] emission and SFR can be easily un¬ 
derstood both in PDRs, where the emission of radiation is proportional to the strength of the far 
UV radiation field thus to the SFR, and in ionized regions since their size is proportional to the 
ionization rate and thus to the SFR. 

For redshifted [CII] line, the redshift range corresponds to observing frequencies of ~ 200 — 
300 GHz. We assume a future [CII] intensity mapping instrument, based on a grating spectrometer 
and 20,000 bolometer array detectors for spectral line measurements. The instrument is assumed 
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to be on a telescope with a 10-m aperture. The details specifications are listed in Table 2. The 
forecasted cross-power spectrum of 21cm and [CII] surveys at z = 8 is shown in Figure 2. The red 
curve is the predicted amplitude of cross-correlation based on [CII] models in Gong et al. (2012), 
while the green curves indicate the theoretical uncertainties. Note the cross-correlation is negative 
at k <~ 1 (h/Mpc), a signature from the anti-correlation of [CII] and 21cm on scales larger than 
the typical bubble size, as the former traces star formation activities thus the ionized region, while 
21cm traces the neutral part of the Intergalactic Medium (IGM). A 4.7-a detection is expected 
with this particular setup. At these frequencies, however, contributions from different CO rota¬ 
tional line emissions coming from different redshifts may confuse the redshifted [CII] emissions. 
One may apply bright-source masking or template fitting techniques to extract the redshifted [CII] 
signals. A pilot [CII] intensity mapping experiment, Time-Pilot, is currently underway to map out 
the redshifted [CII] emissions from high redshifts (Crites et al. 2014). 



k (h/Mpc) 


Figure 2: The cross-power spectmm of SKA1-LOW 21cm with a potential [CII] line mapping program at 
z = 8. The expected signals are plotted in red, while the theoretical uncertainties of the models are indicated 
by the green curves. The error bars are calculated based on parameters listed in Tables 1 and 2. 


4. Lya Intensity Mapping 

Lya photons have a rest-frame wavelength of 1216 A and so Lya emission during the EoR 
will redshift to the near-infrared regime today, making it potentially detectable by narrow-band 
infrared detectors. Lya photons emitted by galaxies are mostly absorbed and reemitted by the 
neutral hydrogen in the galaxy which causes a scatter of the radiation greatly decreasing the Lya 
flux detected by direct observations of Lya emitters. Therefore, galaxy surveys are not able to fully 
measure all of the intrinsic Lya emission. Intensity mapping is however a low resolution technique 
and so by not attempting to resolve the sources of Lya photons we can in principle detect all of the 
Lya radiation emitted both from galaxies and from the IGM. 
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Table 2: Experimental Parameters for a Possible [CII] Mapping Instrument. 


Aperture diameter (m) 

10 

Survey Area (As; deg 2 ) 

16 

Total integration time (hours) 

4000 

Tree spectral range ( B v ; GHz) 

185-310 

Lreq. resolution ( S v ; GHz) 

0.4 

Number of bolometers 

20,000 

Number of spectral channels 

312 

Number of spatial pixels 

64 

Beam size a (Obeand LWHM, arc min) 

0.4 

Beams per survey area 3 

2.6 x 10 5 

<T p j x : Noise per detector sensitivity 2 (Jy^/s/sr) 

2.5 x 10 6 

/p’ 1 ]): Integration time per beam 2 (hours) 

1.0 

z = 8 Vpix (Mpc/h) 3 

4.8 

z = 8 P^ 11 (Jy/sr) 2 (Mpc/h) 3 

4.3 xlO 9 


2 values computed at 238 GHz, corresponding to [CII] at z = 7. 


Recent work by Silva et al. (2013); Pullen et al. (2014) describes the process of Lya emission 
in the EoR and post-EoR and shows estimates for intensity mapping of the Lya signal at redshifts 
7 < z < 11. Lya emission from galaxies is mainly sourced by stellar ionizing radiation since stars 
emit photons which ionize neutral hydrogen, which then emits Lya photons upon recombination, 
and also because the heating of the gas by stellar UV radiation gives rise to e-HI collisions causing 
further Lya emission. During the EoR stellar populations can also source emission in the dense and 
ionized IGM surrounding the galaxies through e-p recombinations and e-HI collisions. In addition, 
diffuse Lya emission is also important. It can originate either from e-p recombinations sourced by 
X-ray radiation or from stellar continuum photon redshifted into the Lya line of the IGM. 

The intensity of the contributions from galaxies and from the IGM is dependent on several 
key astrophysical parameters such as: the star formation efficiency, the stellar spectrum, the escape 
fraction of ionizing photons from galaxies to the IGM, the gas temperature and clumping, the 
minimum mass of Lya halos and the ionization state of the IGM. The relative contribution from 
galaxies and from the IGM to the Lya photons budget is highly dependent on the ionization history 
and on the local heating of the IGM and so it is very difficult to estimate. 

Observational maps of Lya emission will be contaminated by extragalactic continuum emis¬ 
sion and line foregrounds and also by emission from our galaxy. Continuum contamination can 
in principle be removed from intensity maps taking into account the smooth evolution of this ra¬ 
diation with frequency compared to the evolution of the Lya line; however, zodiacal light emitted 
from our galaxy will bring confusion to the observational maps making it possible to only extract 
the Lya power spectra at small scales where zodiacal light is spatially smooth. Loreground lines 
from lower redshifts, namely the 6563 A Ha, the 5007 A [OIII] and the 3727 A [Oil] lines will 
strongly affect Lya observational maps however their contamination can be removed by masking 
the contaminated observational pixels as was shown in Gong et al. (2014). 

Below we illustrate the cross-correlation signature of Lya and 21cm emissions during EoR, 
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Table 3: Experimental Parameters for a Possible Lya Mapping Instrument. 


Aperture diameter (m) 

0.2 

Survey Area (As; deg 2 ) 

13 

Total integration time (hours) 

2900 

Free spectral range (By ; /./ m) 

0.85-1.1 

Freq. resolution (A/5^) 

220 

Number of pixels in 2D array 

72900 

FOV per pointing; deg 2 

0.6 

Observational time per pointing (hours) 

129.5 

Survey volume (Mpc/h) 3 

8.5 x 10 7 


from an assumed Lya intensity mapping experiment which consists of an aperture array with the 
parameters described in table 3. The parameters of the 21cm intensity mapping observation with 
SKA1-LOW are described in Table 1 for z = 8, and the assumed frequency dependence of instru¬ 
ment system temperature T sys and the collecting area A e are as follows: T sys = T^ y + T rec , where 
T s kv = 60 ( 300I ^ Hz j 2 ' 55 K is the sky temperature in Kelvin, and T rec = 0.1* T s ^ y + 40 K is the instru¬ 
ment receiver temperature. A e = 925 (-yr)Nn 2 . The Lya intensity mapping calculation is from 
Silva et al. (2013). The forecast shows that the 21cm and Lya cross-power spectra can be de¬ 
tected at high SNR values of (789, 442, 273, 462) for z = (7,8,9,10), respectively, and is shown in 
Figure 3. 



k[h/Mpc] 

Figure 3: The SKA1-LOW 21cm and Lya cross-power spectra at z = (7,8,9,10). The Lya models are 
based on Silva et al. (2013). The cross-power spectra are predicted to be detectable at > 100 SNR signifi¬ 
cance level, given the assumed survey parameters 
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5. Discussion 

The strength of the correlation between HI and molecular line emission depends on multiple 
factors and in particular, the sign of the correlation will depend on whether the emission comes 
mostly from the galaxies or the IGM. When cross-correlating CO with HI, the cross-correlation 
signals are expected to be associated with the clustering of galaxies (CO) and the IGM (HI), thus 
there is a strong and easy to interpret anti-correlation. When cross-correlating [CII] with HI, the 
bulk of the [CII] emissions are expected to come from galaxies, with some small contributions 
from the IGM, and we expect a strong anti-correlation. At the high redshifts of interest (z>6), a 
small amount of metals may reside in the IGM, however, the [CII] spin temperature is expected 
to follow the CMB temperature (Gong et al. 2012) thus little [CII] emission from the IGM is 
expected. For cross-correlations of HI with Lya emissions, the situation is more complicated since 
the Lya emission from the IGM can be very high and is very uncertain. Depending on the model 
considered, the emission from the IGM can even be higher than the emission from galaxies at some 
redshifts. In this case, Lya and HI would be positively correlated on large scales; on the other hand, 
if no Lya-HI cross-correlation signals are found, we can place constraints on the Lya sources of 
emission. 

6. Summary 

In this chapter we motivated a novel use of SKA1-LOW 21cm to probe the epoch of reion¬ 
ization. By cross-correlating the 21cm line with other atomic and molecular lines observed in the 
intensity mapping regime, we can not only validate a potential 21cm EoR detection but also learn 
more about the EoR itself. On the one hand, the 21cm line traces neutral regions, the yet to be 
ionized universe while the other atomic lines trace star formation activities. The combination thus 
forms a potent, complete and unique picture of the reionization process. 

We focus our study on CO, [CII] and Lya transition lines currently identified as the most 
promising tracers. For CO, we study the possibility of using directly SKA1-MID at appropriate 
frequencies to map the EoR at z>7.5. For [CII] and Lya, we assume the successful deployment of 
relevant instruments currently being planned. In all cases, we find that the strength of our detection 
strongly varies with theoretical model. For example, in the case of CO, according to the particu¬ 
lar models considered, we could go from a strong detection (greater than 6 a) to a non-detection. 
The same holds for [CII], and it appeal's more promising for Lya to be detected at high SNRs . 
While this situation might be worrisome, it simply reflects the fact that we are probing a totally 
new territory for astrophysics. Besides, these line tracers have different sources of astrophysical 
contamination which are in general much less severe than the ones plaguing the redshifted 21 cm 
line. A cross-correlation measurement can thus serve as an independent confirmation of the cos¬ 
mological origin of the measured signals. The use of multiple line tracers would thus be invaluable 
to validate and enrich our understanding of the EoR, and it will open up a huge discovery space. 
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